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Mellin-Bai^nes and sector decomposition methods are used to evaluate tensorial Feynman diagrams in the Euclidean kinematical 
region. Few software packages are shortly described and few examples demonstrate their use. 



1. Introduction 



2. General tensor structure 



Mellin-Barnes (MB) methods (see e.g. fT\ and refer- 
ences therein) and sector decomposition (SD) methods 
(see e.g. f2] and references therein) are frequently 
used in theoretical high energy physics calculations. 
We will restrict ourselves here mainly to Feynman loop 
integrals in d dimensions. There are nowadays public 
software packages which may help to automatize these 
calculations. On the MB Tools webpage 1 3 1 there are 
useful codes which realize: analytic continuations of 
Mellin-Barnes integrals - MB.m |4| and MBresolve.m 
|[5ll ; parametric expansions of Mellin-Barnes integrals - 
expansion. m (author M. Czakon) andbarnesroutines.m 
(author D. Kosower). For the creation of MB repre- 
sentations and later calls of other tools, AMBRE.m 
can be used |6 7 |, which we describe here. AMBRE 
vl.x generates multiloop scalar [planar] and one-loop 
Feynman integrals. Here we extend the package to 
v2.0, which may treat also tensor [planar] structures. 
For the sector decomposition method, there are two 
public programs. FIESTA 1 8 1 and FIESTA2 1 9 1 calcu- 
late scalar Feynman integrals with Mathematica, while 
Bogner/Weinzierl's (BW) program 111 01 allows to cal- 
culate basic polynomials needed for the calculation of 
any Feynman integral with the C++ based computer al- 
gebra system Ginac j 1 1 j. Here we describe the Mathe- 
matica program CSectors. m which may serve as a kind 
of user interface to BW and Ginac in order to calculate 
fully automatically tensor Feynman integrals. 



*Talk presented by J.G. at "Loops and Legs in Quantum Field The- 
ory", 10th DESY Workshop on Elementary Particle Theory, 25-30 
April 2010, Worlitz, Germany 



We consider the L loop Feynman momentum inte- 
grals with propagators P„ and tensor structure T{k): 
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The momenta qt are linear combinations of internal 
momenta kj and external momenta p„, and T{k) is 
a tensor structure, e.g. {\^k^ jk^k'^^,. . .). Further, 
A^v = L ^n- The general scalar L loop integrals Gl[1] 
are implemented in AMBRE.m vl.x, however tensor 
structures were allowed only for the one-loop case, 
G\[T{k)]. The symbol {R}V'v'h] means the totally 
symmetrized sum of all possible combinations in the 
argument. We give as an explicit example the rank 
R = 3 tensor: 

_ pMl pM2 pM3 _|_ ^Mi M2 pM3 _|_ ^M2M3 pMi 
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It is Aq — — 1, and A,- is zero for r odd, 
while = for r even. Further, 

= and p!^' = LiiMaiQi]^'- The 

M and Qi are defined from ELiPiXi = - 

= fefM,y^,- - 2I^^i k^jQj + 7 and M = 

det(M)M^'. The C/ and F polynomials were defined 
e.g. in J6l7l . 

3. The package AMBRE.m 

In 1121 it is shown how the tensor structures intro- 
duced in section|2]lead, in the loop-by-loop approach, 
to a list of Mellin-Barnes representations. For that pur- 
pose, the following basic function of the AMBRE.m 
package is invoked: 

MBrepr[{numerator},{propagators}, 
{internal momenta}]. 

Let us take the example of figure [T] and sketch a 
sample math file. 

First, packages must be loaded: 

« MBvl.2.m 
« AMBRE.m 

« MBnum.m (3) 

or 

« MBresolve.m (4) 
« barnesroutines .m 

As indicated, barnesroutines .m works with ver- 
sion 1.2 of MB.m. If we want to use this package, the 
MB integrals up to the needed power in e must be pre- 
pared first (by analytic continuation), and the second 
part of the sample MB_SE510m.m file may look as 
follows: 

invariants = {pl~2 -> s}; 

repr = MBrepr [{kl*pl , kl*pl, k2*pl}, 
{PR [kl , , nl] *PR [k2 , , n2] *PR [k2-kl , , n3] 
*PR [kl+pl , , n4] *PR [k2+pl , , n5] } , 
{k2,kl}] ; 

in the case of the MBnum.m package, followed by: 

^Kinematics in tiie form of a list must be also properly defined. 
"^See also the complete MB_SE510m.m source at [7 1. 



SetOptions [MBnum, Analytical -> True, 
Numerical -> False, ShowMBrep -> True]; 

MBanalytic = MBnum [repr, 1, {s -> -11}, 
{nl->l, n2->l, n3->l, n4->l, n5->l}, 2]; 

(5) 

or, in the case of the MBresolve.m package, fol- 
lowed by: 

stepl=MBresolve [#/ . powers , eps] &/ Srepr ; 
step2=MBexpand [stepl , 

Exp [2*eps*EulerGainma] , {eps, 0,1}] ; 
MBanalytic=MBmerge [step2] ; 

(6) 

Finally, after getting the MB integrals at the e level 
needed (here second argument of the MBnum func- 
tion), we "Simplify" the result using Barnes lemmas 
with the internal MB.m functions MBmerge and MBin- 
tegrate: 

after = 

Process [MBanalytic , Range [10] ] //MBmerge ; 
SEntimMB=MB integrate [after ,{s->-ll}] 

The phrase Rcinge [10] means that we assume to get 
up to 10 dimensional MB representations; in fact, they 
are maximally 2-dimensional here. 
The final result is: 

SEnumMB = -7.5625/e2 - 20.4506/e 
-178.18±0.0171936 
+(18.3642 ±0.0248465)e. (7) 

Singularities are here at most 1 -dimensional integrals, 
so MB.m solves them within Mathematica very accu- 
rately and no error is given. 

4. The package CSectors.m 

The package CSectors.m allows to derive integral 
representations for Feynman integrals, which are de- 
composed into so-called sectors with a well-defined 
singularity structure. The invariants should be of the 
Euclidean type. The structure of the main function is 
analogous to that in Ambre.m: 

^Actually the use of barnesroutines. m does not help much here, and 
it takes most of the running time: the time of the total calculation is 
17 seconds without applying Barnes Lemmas. 
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Figure 1 . A 2-loop self energy topology. 



DoSectors[{numerator},{propagators}, 
{internal momenta}] [low,des] 

The sample executable file might look as follows: |^ 

# ! /bin/bash 
math « \FILE 

<< CSectors.m 

invariants = {pl"2->s} 

SetOptions [DoSectors , SetStrategy -> C, 
TextDisplay->True , SourceNaine->SE510m_num, 
TempFileDelete->False] ; 

max=0; s=-ll; 

invariants = {pl~2 -> s}; 

SEnumSD = DoSectors [{kl*pl ,kl*pl ,k2*pl}, 
{PR [kl , , 1] *PR [k2 , , 1] *PR [k2-kl ,0,1]* 
PR [kl+pl ,0,1] *PR [k2+pl , , 1] > , {k2 , kl>] 
[-2 ,max] 

FILE 

The numerical result in the example is (see also file 
SD_LL2010.out at 1131 ): 

SEnumSD = -178.1927-7.56258/e^ 
-20.4505/e + 18. 3 94000000000005 e, 
{0.01 1 130644628528934, 0.00029563 /e^, 
0.00302109/e,0.06629380324170578*e}. (8) 

The second bracket gives an error estimation. Due 
to the tensor structure (|T|l, a number A of integrals /,■ 

■•See file SD_SE510m.sh at (13). 




Figure 2. The planar 2-loop Bhabha topology Bl. 



appear at a given e level. An error for the numerical 
result Ig at £ level is estimated in a standard way: 




In the MB.m package the error AI is, for the single MB 
integral /, controlled by switching on an appropriate 
MB.m option: Debug->True. Let us just note that the 
error can be underestimated, see i4j and MB.m for 
details. 

5. Mellin-Bames approach versus Sector Decom- 
position approach 

All the results given in the last section in ([7]i and ([8]l 
have been obtained with default parameters needed for 
numeric defined by the WB and MB.m packages. The 
total time of calculation on a Xeon PC with 32 GB 
RAM was in the MB case 17s and about 5 minutes for 
both the C- and X-strategies in the SD case. 

Calculational times can differ substantially depend- 
ing on methods and integrals. In table [T] and table [2] 
results for the massless and massive planar 2-loop 
topology Bl are given, see Fig.[2]and 1 14 15 1 for the 
topology's definition. 

We can see that in the case of the massless box Bl, 
the MB method can be much faster than the SD method. 
In the massive case the difference is smaller, moreover, 
the SD method can be often faster. In addition results 
obtained by FIESTA2 |9 1 are presented. Note a differ- 
ence in sign between FIESTA2 and MB/SD, see JSJ 
for FIESTA2 conventions. 

For completeness, the FIESTA2 sample Mathemat- 
ica file for getting the numbers in Table |2] is repro- 
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Table 1 

The scalar massless 2-loop planar double box Bl, 
shown in figure |2] The scalar case allows for direct 
comparisons of the packages discussed here with FI- 
ESTA2. For source and output files see [El and M13I . 



massless 


AMBRE and MB 


e" 


-0.1034 ±6- 10"'^ 


£-1 


-0.10907 


£-2 


-0.00966 


£-3 


0.083188 


£-4 


-0.022857 


T [s] 


22 


massless 


CSectors, X-strat. 


£'' 


-0.1035 ±0.0002 


£-1 


-0.10915±0.00008 


£-2 


-0.00966 ±0.00001 


£-3 


0.083191 ±0.000005 


£-4 


-0.0228574±M0-'' 


T [s] 


1712 


massless 


FIESTA2 


£" 


0.103267 ±0.001 117 


£-1 


0.109049 ±0.000348 


£-2 


0.009656 ±0.000082 


£-3 


-0.08319 ±0.000012 


£-4 


0.022858 


T [s] 


350 


i= -5, f = -7 



duced: 

« FIESTA_2.0.0.m; 
CIntegratePath=" . /CIntegrateMP" 
QLinkPath = "./qLink64" 

invariants = 

{pl'2->m~2 , p2'2->m'~2 , p3'2->m~2 , p4'2->m~2 , 
pl*p2->l/2*s-in~2,p3*p4->l/2*s-m'2, 
pl*p3->l/2*t-m~2 , p2*p4->l/2*t-m'2 , 
p2*p3->l/2*u-m'~2,pl*p4->l/2*u-m~2> 
/.u->4*m'2-s-t/.m->l//Simplify; 

kinematics = {s -> -5, t -> -7}; 

integral=PR[kl , 1 ,nl] *PR [kl+pl ,0,n2] * 
PR [kl+pl+p2 , 1 , n3] *PR [kl-k2 , , n4] * 
PR [k2 , 1 , n5] *PR [k2+pl+p2 , 1 , n6] * 
PR[k2 + pi + p2 + p4, 0, n7] ; 



Table 2 

The scalar massive 2-loop planar double box Bl. 



massive 


AMBRE and MB 


£" 


0.2246 


£-1 


0.06359 


£-2 


-0.023524 


T [s] 


42 


massive 


CSectors, C-strat. 




0.2246 ±0.0001 


£-1 


0.06357 ±0.00003 


£-2 


-0.023524±4-10-^ 


T [s] 


362 


massive 


FIESTA2 


£" 


-0.224756 ±0.000485 


£-1 


-0.063622 ±0.000102 


£-2 


0.023528 ±0.00001 


T [s] 


35 


s = 


—5, / = —7, m = 1 



f integral=(List @@ integral) 
/.PR[q_,m_,n_]->-(q~2-m~2) ; 

SDEvaluate [UF [{kl , k2} , f integral , 

Join [invariants , kinematics]], 
{1,1,1,1,1,1,1},0] 

6. Summary 

The CSectors. m package has been prepared in 
order to allow an easy use of the SD package by 
Bogner/Weinzierl for the evaluation of Feynman ten- 
sor integrals. The AMBRE package has been enlarged 
to involve tensor structures for multiloop integrals in 
an automatic way. We are planning for the near future 
to automatize the construction of MB representations 
for nonplanar diagrams, what deserves not to use the 
loop-by loop approach. Further, we will foresee to 
build MB representations for special forms of linear 
propagators which are present e.g. in the HQET theory 
or in calculations of the QCD static potential (see the 
talk by V. Smirnov at this conference 0161 ). 
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